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TECHNICAL  NOTE  5179 

A THEORETICAL  INVESTIGATION  OF  THE  HEATING-UP  PERIOD 
OF  INJECTED  FUEL  DROPLETS  VAPORIZING  IN  AIR 
By  M.  M.  El  Wakil,  0.  A.  Uyehara,  and  P.  S.  Myers 


SUMMARY 


Theoretical  calculations  were  performed  for  the  primary  purpose  of 
determining  whether  the  unsteady  state,  that  is,  the  time  during  which 
the  droplet  is  heating  up,  or  cooling  down,  to  the  wet-bulb  temperature, 
is  of  importance  when  single  fuel  droplets  are  vaporizing  under  condi- 
tions of  interest  in  jet  engines.  The  results  are  presented  as 
temperature -time,  mass-time,  and  penetration-time  histories  as  well  as 
heat  and  mass  transfer  histories  of  single  droplets.  These  calculations 
were  conducted  in  the  unsteady  or  nonequilibrium  stage  beginning  at  the 
time  the  single  droplets  come  into  contact  with  the  air  and  up  to  the 
point  at  which  they  reach  the  equilibrium  or  wet-bulb  temperatures,  at 
which  time  all  the  heat  transferred  from  the  air  to  the  droplet  is  carried 
away  by  the  diffusing  vapor.  After  this  point  is  reached,  the  heat  and 
mass  transferred  to  and  away  from  the  droplet  become  functions  only  of 
the  radius  and  velocity  of  the  droplet. 

Calculations  were  performed  for  droplets  of  different  size  and  fuel 
composition  under  different  ambient  air  conditions.  Representative 
results  are  presented  in  this  report.  A complete  set  of  graphs  presenting 
the  results  of  all  the  calculations  is  contained  in  the  manuscript  copy  of 
this  report.  This  copy  is  available  for  loan  or  reference  in  the  Division 
of  Research  Information,  National  Advisory  Committee  for  Aeronautics, 
Washington,  D.  C. 

The  calculations  show  the  unsteady-state  portion  of  the  total  vapori- 
zation time  to  be  of  importance  in  most  cases.  At  high  air  temperatures 
and  for  high-volatility  fuels,  the  unsteady  state  represents  the  greater 
part  of  the  vaporization  time.  The  larger  droplets  would  be  expected  to 
reach  the  combustion  zone  while  still  in  the  unsteady  state. 

In  most  cases,  during  the  unsteady  state,  more  of  the  heat  arriving 
at  the  surface  of  the  liquid  droplet  goes  to  sensible  heat  than  goes  to 
latent  heat  of  vaporization.  This  is  especially  true  for  low-volatility 
fuels  where  the  net  heat  transfer  from  the  air  to  the  film  goes  mostly 
into  heating  up  the  liquid.  The  rate  of  mass  transfer  usually  increases 
with  time  as  the  droplet  heats  up.  At  high  air  densities,  however,  the 
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mass  transfer  decreases  with  time,  although  the  droplet  is  being  heated 
up.  The  time  required  by  a certain  size  droplet  to  reach  its  wet-bulb 
temperature  is  less  for  higher  volatility  fuels  at  the  same  ambient  and 
initial  conditions. 

It  should  be  realized  that  the  calculated  curves  depend  upon  the 
use  of  heat  and  mass  transfer  correlations  that  are  extended  far  beyond 
the  conditions  under  which  they  were  obtained.  Thus,  extensive  use  of 
the  values  shown  on  the  calculated  curves  should  not  be  made  until  either 
the  curves  or  the  extrapolations  of  the  heat  and  mass  transfer  correla- 
tions have  been  checked  experimentally. 

In  addition  to  the  theoretical  work,  in  order  to  verify  the  results 
of  the  calculations,  which  showed  the  unsteady  state  to  be  of  importance, 
apparatus  was  set  up  to  obtain  experimental  temperature -time  histories 
of  comparatively  large  droplets  vaporizing  in  air  streams  at  constant 
velocities  and  atmospheric  pressure.  This  apparatus  is  described  and  a 
few  experimental  curves  are  included.  These  curves  show  the  relative 
distribution  of  the  time  interval  between  the  unsteady  and  steady  states 
when  different  air  temperatures  and  different  fuels  are  used. 


INTRODUCTION 


Knowledge  of  the  phases  through  which  a heterogeneous  mixture  of 
air  and  fuel  passes  before  ignition  takes  place  is  essential  for  the 
rational  design  of  jet-engine  combustion  chambers.  For  discussion  pur- 
poses, the  time  elapsed  before  ignition  is  often  divided  into  physical 
and  chemical  components,  although  the  relative  magnitudes  of  these  two 
components  are  not  known.  Since  the  physical  component  inevitably  pre- 
cedes and  perhaps  in  part  overlaps  the  chemical  component,  knowledge  of 
the  magnitude  of  and  the  factors  affecting  the  physical  component  is  of 
obvious  importance  in  any  study  of  the  time  interval  before  ignition. 

Fuel  injected  through  a nozzle  into  a combustion  chamber  is  believed 
to  pass  through  the  following  stages:  The  fuel  leaves  the  nozzle  orifices 

as  a ligament  or  sheet.  This  ligament  or  sheet  then  breaks  down  into  dif- 
ferent size  droplets.  Appreciable  heat  transfer  from  the  air  present  in 
the  combustion  chamber  to  the  fuel  occurs  only  after  the  ligament  has 
broken  up  into  droplets.  This  is  due  to  the  appreciable  increase  of  sur- 
face area  of  the  fuel  that  is  exposed  to  the  air.  As  heat  transfer  takes 
place,  the  droplets  heat  up  (or  possibly  cool  down,  depending  upon  ambient 
air  conditions  and  initial  fuel  temperature)  and  at  the  same  time  lose 
part  of  their  mass  by  vaporization  and  diffusion  into  the  air,  with  both 
heat  and  mass  transfer  being  markedly  affected  by  droplet  size  and  velocity 
relative  to  the  air.  At  the  same  time,  the  droplets  are  being  slowed  down 
(or  speeded  up,  depending  upon  the  magnitude  of  the  air  and  fuel  velocities) 
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relative  to  the  air  by  aerodynamic  drag  forces.  After  a certain  time 
has  elapsed,  each  droplet  attains  an  asymptotic  or  equilibrium  tempera- 
ture equal  to  the  wet-bulb  temperature  corresponding  to  the  conditions 
present  at  that  moment. 

The  larger  droplets  are  slower  in  attaining  equilibrium  conditions; 
but,  although  they  have  essentially  the  same  initial  speed  as  the  smaller 
droplets,  they  are  slowed  down  relative  to  the  air  at  a lesser  rate,  and 
they  move  along  with  the  air  ahead  of  the  smaller  ones  that  were  injected 
at  the  same  instant.  The  smaller  drops,  however,  give  away  their  mass 
and  completely  vaporize  faster  and  travel  a shorter  distance  through  the 
combustion  chamber  than  the  larger  droplets. 

A cloud  of  vapor  due  to  the  smaller  droplets  is  thus  rapidly  formed 
and  moves  along  with  the  air.  The  mass  of  vapor  given  away  by  the 
incoming  larger  droplets  is  added  to  this  vapor  cloud.  Somewhere  in  the 
combustion  chamber  a combustible  mixture  of  air  and  fuel  vapor  is  formed 
and  is  ready  for  ignition  at  that  point.  In  the  absence  of  outside 
ignition,  this  mixture  must  be  heated  by  the  air  until  it  reaches  its  self- 
ignition temperature  if  combustion  is  to  take  place.  By  definition,  the 
physical  portion  of  the  ignition  delay  extends  until  a combustible  mixture 
is  formed  and  heated  to  the  ignition  temperature.  Since  some  chemical 
reaction  may  take  place  at  temperatures  lower  than  the  ignition  tempera- 
ture, the  chemical  portion  of  the  ignition  delay  period  probably  overlaps 
the  physical  portion.  The  vaporization  process  is  thus  an  essential  part 
of  the  ignition  delay  period. 

Although  the  fuel  droplets  reach  their  equilibrium  or  wet-bulb  tem- 
peratures asymptotically  with  time,  the  vaporization  process  can  roughly 
be  separated  into  the  unsteady  state  and  the  steady  or  equilibrium  state. 
Numerous  investigations  have  been  conducted  covering  the  latter  stage 
(refs.  1 and  2),  but  only  a few  attempts  have  been  made  to  study  the 
unsteady  state  (ref.  3) • It  was  the  intent  of  this  investigation  to 
determine  theoretically  temperature-time  histories  of  droplets  of  dif- 
ferent size  and  fuel  composition  under  different  ambient  air  conditions 
in  an  attempt  to  determine  whether  the  unsteady-state  period  is  of 
importance  in  the  history  of  the  fuel  droplets  previous  to  ignition. 

This  report,  then,  presents  a detailed  study  of  the  unsteady-state 
portion  of  the  total  vaporization  time  of  single  fuel  droplets  injected 
into  air.  Calculations  were  performed  for  numerous  combinations  of  the 
following  conditions: 

Fuels  n-octane,  benzene,  n-pentane 

Initial  radii,  p 10,  50,  100 

Initial  droplet  temperatures,  °F 50,  200 

Initial  droplet  velocities,  fps  .....  100,  400 

Air  pressures,  atm  0.5>  1,  5 

Air  temperatures,  °F 100,  500,  1,000 


If 
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The  three  fuels  have  standard  atmospheric  toiling  temperatures  of  96.9°, 
176.2°,  and  258.2°  F,  respectively.  The  calculations  covered  all  com- 
binations of  the  above  conditions  except  when  the  initial  droplet  tem- 
peratures were  higher  than  the  boiling  points  of  the  fuels  at  the  corre- 
sponding air  pressures . 

The  properties  of  the  three  fuels  under  investigation  and  for  air 
were  obtained  from  various  references  as  shown  in  table  I. 

This  work  was  conducted  at  the  University  of  Wisconsin  under  the 
sponsorship  and  with  the  financial  assistance  of  the  National  Advisory 
Committee  for  Aeronautics. 


SYMBOLS 


Ac  cross-sectional  area  of  liquid  droplet,  sq  in. 

A0  surface  area  of  liquid  droplet,  sq  in. 

B0  thickness  of  air-vapor  film  surrounding  droplet,  in. 

Cp  coefficient  of  drag  for  spheres,  unitless 

Cp  specific  heat  of  air  at  constant  pressure,  Btu/(lb)(°F) 

*Ei 

Cr,^  specific  heat  of  fuel  vapor  at  constant  pressure, 

Pf  Btu/(lb)(°F) 

Cp^  specific  heat  of  liquid  fuel,  Btu/(lb)(°F) 

Cp  average  specific  heat  of  air-vapor  mixture  in  film  at 

m constant  pressure,  Btu/(lb)(°F) 

Dv  diffusion  coefficient  of  air-vapor  system,  sq  in. /sec 

d penetration  or  distance  traveled  relative  to  air,  in. 

F drag  force,  lb 

ga  weight  fraction  of  air  in  film,  unitless 

gp  weight  fraction  of  fuel  vapor  in  film,  unitless 

h coefficient  of  heat  transfer  in  film  in  absence  of  mass 

transfer,  based  on  surface  area  of  liquid  droplet  Aq, 

Btu/(sq  in. ) (sec) (°F) 
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hL 

Ki>  k2, 
ka 

kf 

kg 

kL 

M 

Ma 

% 


M, 


lm 


m 


N, 


Nu 

%u' 


■Tr 


■vRe 


N. 


Sc 


PaB 


coefficient  of  heat  transfer  to  liquid,  Btu/(sq  in.) (sec) (°F) 
Kj,  . . . constants 

thermal  conductivity  of  air,  Btu/(in.)(sec)(°F) 

thermal  conductivity  of  fuel  vapor,  Btu/(in. ) (sec) (°F) 

coefficient  of  mass  transfer  in  film,  l/sec 

thermal  conductivity  of  liquid  fuel,  Btu/(in. ) (sec) (°F) 

average  thermal  conductivity  of  air-vapor  mixture  in  film, 
Btu/( in. ) (sec) (°F) 

ratio  of  instantaneous  mass  to  original  mass  of  droplet, 
percent 

molecular  weight  of  air,  lb/mole 

molecular  weight  of  fuel  vapor,  lb/mole 

apparent  molecular  weight  of  air-vapor  mixture,  lb/mole 

mass  of  liquid  droplet,  lb 

Nusselt  number  for  heat  transfer,  unitless 

Nusselt  number  for  mass  transfer,  unitless 

Prandtl  number,  CpmUrn^m>  unitless 

Reynolds  number,  2r0Up^i%,  unitless 


Schmidt  number 


Pm/pmDv> 


unitless 


total  pressure,  lb/sq  in.  abs 

average  partial  pressure  of  air  in  film,  Pj,  - 
lb/sq  in.  abs 


PfT 


partial  pressure  of  air  at  outer  edge  of  film, 
lb/sq  in.  abs 


partial  pressure  of  air  at  liquid  droplet  surface, 
lb/sq  in.  abs 

average  partial  pressure  of  fuel  vapor  in  film, 
lb/sq  in.  abs 

partial  pressure  of  fuel  vapor  at  outer  edge  of  film, 
considered  to  be  zero,  lb/sq  in.  abs 

partial  pressure  of  fuel  vapor  at  liquid  droplet  surface, 
equal  to  vapor  pressure  of  fuel  at  droplet  tempera- 
ture T^,  lb/sq  in.  abs 

total  heat  transfer  from  air  to  droplet,  Btu/sec 

sensible  heat  received  by  droplet,  Btu/sec 

heat  carried  away  from  droplet  by  diffusing  vapor  in  form 
of  superheat,  Btu/sec 

heat  received  at  droplet  surface,  Btu/sec 

heat  of  vaporization,  Btu/sec 

universal  gas  constant,  in-lb/(mole) (°F) 

radius  at  any  point  in  film,  in. 

molal  rate  of  diffusion  of  air  in  equimolal  diffusion, 
moles/(sec) (sq  in.) 

molal  rate  of  diffusion  of  fuel  vapor  in  equimolal 

diffusion,  based  on  surface  area  of  liquid  droplet  Aq, 
moles/(sec) (sq  in.) 

radius  of  liquid  droplet,  in. 

temperature  in  film  at  radius  r,  °R 

asymptotic  or  wet-bulb  temperature  of  droplet,  °R 

temperature  of  air  outside  film,  °R 
temperature  of  liquid  droplet,  °R 
logarithmic  mean  temperature  in  film,  °R 
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U droplet  velocity  relative  to  air,  in. /sec 

w mass  of  fuel  vapor  diffused  out,  lb/sec 

x proportionality  factor 

r - ro 

y fraction  of  film  thickness  at  radius  r,  , unitless 

Bo 

z factor,  wCp^hA0,  unit  less 

a correction  factor  of  mass  transfer  due  to  semipermeability 

of  film,  unities s 

0 time,  sec 

A latent  heat  of  vaporization  of  fuel,  Btu/lb 

pa  viscosity  of  air,  lb/(in.)(sec) 

Pf  viscosity  of  fuel  vapor,  lb/(in. ) (sec) 

pm  average  viscosity  of  air-vapor  mixture  in  film,  lb/(in. ) (sec) 

Pa  density  of  air,  lb/cu  in. 

0^  density  of  liquid  droplet,  lb/cu  in. 

Pm  average  density  of  air-vapor  mixture  in  film,  lb/cu  in. 

ANALYSIS 


In  the  following  theoretical  study  of  the  vaporization  of  single 
droplets,  two  assumptions  were  made.  The  first  assumption  was  that  an 
endless  supply  of  air  surrounded  the  droplet  at  all  times.  This  means 
that  there  is  no  interaction  between  different  droplets  in  the  same  air 
stream,  or  that  only  a single  droplet  was  considered  present.  The 
partial  pressure  due  to  the  fuel  vapor  at  the  outer  edge  of  the  air- 
vapor  film  surrounding  the  droplet  was  thus  considered  to  be  zero  at  all 
times.  The  second  assumption  was  that  the  air  pressure  and  temperature 
remained  unchanged  throughout  the  vaporization  process.  It  is  known 
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that  during  injection  in  a continuous-flow  air  stream  the  air  temperature 
decreases  because  of  heat  transfer  and  that  in  the  absence  of  friction 
its  pressure  should  increase  because  of  the  slowing  down  of  the  stream; 
the  second  assumption  is  consistent,  however,  with  the  assumption  of  a 
single  droplet  in  a large  volume  of  air. 

Figure  1 shows  a spherically  symmetrical  droplet  surrounded  by  a 
film  containing  an  air  - fuel- vapor  mixture.  At  any  one  instant,  the 
mixture  strength  varies  from  a maximum  at  the  droplet  surface  to  a mini- 
mum at  the  outer  edge  of  the  film.  The  shape  of  the  temperature  gradient 
curve  depends  on  the  relative  temperatures  of  the  fuel  droplet  and  air 
and  the  mass  of  vapor  being  diffused  out  at  that  instant.  Under  the 
ambient  conditions  considered  in  this  study,  eliminating  the  case  where 
the  droplet  was  introduced  with  a vapor  pressure  higher  than  the  total 
pressure,  the  vaporization  process  can  take  place  in  three  different 
cases: 

(1)  When  the  initial  temperature  of  the  droplet  is  lower  than  the 
air  temperature  and  lower  than  its  own  wet-bulb  temperature,  in  which 
case  the  droplet  heats  up  continuously  during  the  unsteady-state  portion 
of  the  vaporization  time.  This  case  comprises  the  larger  part  of  the 
calculated  conditions  contained  in  this  report. 

(2)  When  the  initial  temperature  of  the  droplet  is  lower  than  the 
air  temperature  but  higher  than  its  own  wet-bulb  temperature  at  the 
ambient  air  conditions.  In  this  case  the  droplet  cools  down  continu- 
ously while  vaporizing  and  the  sensible  heat  lost  by  the  liquid  supplies 
part  of  the  heat  of  vaporization  of  the  diffusing  vapor.  Heat  is  always 
transferred  from  the  air  to  the  droplet,  however. 

(3)  When  the  initial  temperature  of  the  droplet  is  higher  than  both 
the  air  temperature  and  its  own  wet-bulb  temperature.  In  this  case  the 
droplet  cools  down  continuously  until  it  reaches  its  wet-bulb  temperature 
which  is  lower  than  the  air  temperature.  During  that  time  when  the  droplet 
temperature  is  higher  than  the  air  temperature,  there  is  a net  heat  trans- 
fer from  the  droplet  to  the  air  and  thus  heat  and  mass  transfer  take  place 
in  the  same  direction.  When  the  droplet  temperature  has  dropped  below 

the  air  temperature,  the  process  is  the  same  as  that  in  case  (2). 

It  is  convenient  to  divide  the  analysis  into  three  parts  - the  heat 
transfer,  the  mass  transfer,  and  the  velocity  change. 


Heat  Transfer 

Referring  to  figure  1,  and  assuming  that  the  droplet  is  vaporizing 
under  the  conditions  of  case  (l),  the  total  or  net  heat  transfer  Q from 
the  air  to  the  film  surrounding  the  droplet,  prior  to  ignition,  goes  three 
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ways:  (l)  To  heat  up  the  liquid  droplet,  (2)  to  vaporize  the 

liquid,  and  (3)  to  he  carried  or  peeled  hack  with  the  diffusing 

vapor  in  the  form  of  superheat,  Qs.  The  heat  that  arrives  at  the  drop- 
let surface  is  denoted  hy  Qy  and  equals  the  sum  of  and 

It  should  he  noted  that  in  developing  the  relationship  used  Q does 
not  include  radiant  heat  transfer. 

If  at  any  instant  during  vaporization  a point  in  the  air -vapor  film 
at  a radius  r from  the  center  of  the  droplet  and  at  temperature  T is 
considered: 


Qv  = Ql  + % 

= Q - Qs 


(1) 


\ - Mo  f - "Si1  - Tl) 


(2) 


where  w is  the  rate  at  which  the  vapor  is  diffusing  out,  y is  the 
fraction  of  the  film  thickness  at  radius  rQ,  h is  the  coefficient  of 

heat  transfer  through  the  film  in  the  absence  of  mass  transfer,  Cp^  is 

the  specific  heat  of  the  fuel  vapor,  T is  the  temperature  of  the  mix- 
ture at  the  position  y,  and  T^  is  the  temperature  of  the  surface  of 

the  droplet. 

Rearranging  equation  (2)  gives: 


fly 


hA„ 


<*v  “ "S/L  + ^p ? 


dT 


Pf 


Integrating  across  the  film  between  y = 0 to  y = 1.0  where  T = Tj, 
to  T = Tjj  and  considering  the  specific  heat  of  the  fuel  vapor  Cp_^ 
to  he  constant  throughout  the  film: 


hA0  ^v  + ^p-c.  (TB  " tl) 

1 = lo§e  


wCp 


<*v 


\ 
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or  by  rearranging  and  expressing  as  an  exponential: 


jCpf/hA0  Qy  + vcPf (Tb  - Tl) 

Qy 


Defining  z as  wCp^hA0  gives: 


ez  - 1 


Qy  " *%(%  “ \) 


ez  - 1 


(3) 


Qv  = Q 


,eZ  - 1 


(4) 


Thus  the  factor  z /(ez  - 1)  represents  that  fraction  of  the  total 
heat  transfer  Q from  the  air  that  finally  arrives  at  the  surface  of 
the  liquid  droplet  and  supplies  the  latent  heat  of  vaporization  for  the 
vapor  diffusing  out  as  well  as  the  sensible  heat  added  to  the  liquid 

droplet  itself.  The  factor  z/(ez  - l)  therefore  represents  a correc- 
tion factor  to  the  heat  transfer  coefficient  h without  mass  transfer 
from  the  droplet. 

Several  investigators  have  presented  correlations  for  the  heat  trans- 
fer coefficient  h.  However,  as  will  appear  later,  it  was  necessary  for 
these  studies  that  a mass  transfer  coefficient  also  be  available.  The 
correlations  used  in  this  paper  were  taken  from  Ranz  and  Marshall's  work 
(ref.  l).  It  was  felt  that  the  nature  of  the  investigation  conducted  in 
their  work  was  similar  to  the  problems  being  studied  in  this  paper, 
although  the  maximum  air  temperatures  used  in  their  work  were  not  so  high 
as  the  maximum  air  temperatures  used  in  this  investigation. 

From  their  correlation  the  Nusselt  number  for  heat  transfer  used  is: 


Ni 


'Nu  “ 


h(2rp) 


= 2 + 


0.6(Hpr)l/3(%e)l/a 


\ 


(5) 
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The  factors  in  this  equation  were  evaluated  for  the  average  film  con- 
ditions. This  required  knowledge  of  the  air  and  droplet  temperatures 
as  well  as  the  partial  pressures  in  the  film  (see  appendixes  A and  B) . 
The  Ranz  and  Marshall  correlations  were  obtained  from  experiments  con- 
ducted at  low  temperatures  and  low  vapor  pressures.  In  their  work, 

therefore,  the  correction  factor  z /(ez  - l)  was  unity  and  it  is  not 
clear  whether  their  correlation  includes  this  factor  as  a variable. 

Since  the  correlation  was  to  be  used  under  conditions  where  this  factor 
was  not  unity,  and  since  the  correlation  is  similar  in  form  to  those 
obtained  without  mass  transfer,  it  was  decided  to  extrapolate  the  corre- 
lation by  assuming  that  it  did  not  include  the  correction  factor.  The 
correction  factor  therefore  was  computed  and  used  in  all  of  the  calcu- 
lations presented  herein. 

It  can  be  seen  from  the  above  equations  that  a calculation  of  the 
heat  transfer  through  the  film  at  any  instant  during  the  vaporization 
process  of  the  droplet  would  require  the  knowledge  of  the  temperatures 
on  both  sides  of  the  air -vapor  film  surrounding  the  droplet,  of  the 
velocity  of  the  droplet  relative  to  the  air,  and  of  its  radius,  which 
is  dependent  on  the  mass  transferred  up  to  that  point  and  on  the  droplet 
temperature . 


Mass  Transfer 

Neglecting  thermal  diffusion  and  assuming  all  diffusion  to  result 
from  the  driving  force  of  a concentration  or  partial  pressure  gradient 
which  exists  in  the  direction  of  diffusion,  the  following  equations  are 
obtained  for  counterdiffusion  systems  (ref.  4): 


dpf 

dr 


(6) 


and 


iv  =Mf 

xPT 


(7) 


In  this  case  equimolal  rates  of  diffusion  must  take  place  if  the 
total  pressure  is  to  be  maintained,  or: 


rf  = -ra 


(8) 
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and 


PT  = 


pa  + pf 


(9) 


However,  in  most  cases  of  a fuel  droplet  vaporizing  in  air  essen- 
tially unidirectional  diffusion  takes  place.  Equimolal  diffusion  would 
be  approached  only  if  a droplet  approached  its  critical  temperature 
during  the  heating-up  period.  In  the  cases  under  consideration  in  this 
report  essentially  unidirectional  diffusion  exists;  that  is,  there  will 
be  diffusion  of  the  fuel  vapor  away  from  the  droplet  while  the  rate  of 
diffusion  of  the  air  in  the  opposite  direction,  toward  the  droplet,  is 
zero.  In  this  case, 

ra  = 0 (10) 

Combining  equations  (6),  (7),  and  (lO)  and  rearranging: 


dPf  RT 
dr  rfPa 


(11) 


Integrating  between  r = rQ  and  r = rQ  + BQ  where  pa  = PaL 
and  pa  = Pag*. 


r = Dv3?T 
f RTB0 


!oge 


paB 

PaL 


Py 

RTB0 


(12) 


This  is  the  equation  for  the  rate  of  diffusion  of  the  fuel  vapor  in  terms 
of  the  partial  pressures  of  the  air.  Using  equation  (9)  and  assuming  that 
the  partial  pressure  of  the  fuel  vapor  at  the  outer  edge  of  the  film  p^g 

is  zero,  an  assumption  discussed  earlier  in  the  analysis,  the  rate  of  dif- 
fusion of  the  fuel  vapor  can  be  written  in  terms  of  partial  pressures  in 
the  film  as  follows: 
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Dy  PT 

— PfT  ————— 

RTB0  PfL 


logg 


(13) 


where 


rf  = 


RTBr 


PfLa 


a = 


log. 


rT 

PfL 

PT 

PT  " PfL 


(14) 


Equation  (l4)  states  that  the  unidirectional  diffusion  of  the  fuel 
vapor  into  the  air  equals  the  equimolal  rate  of  diffusion,  which  is  a 
function  of  the  concentration,  times  the  factor  a which  is  greater  than 
unity.  This  increase  in  rate  represents  a correction  for  the  fact  that 
the  net  rate  of  movement  of  the  air  is  zero,  while  the  diffusion  veloci- 
ties of  the  fuel  vapor  and  the  air  relative  to  each  other  must  be 
maintained . 

Because  of  the  difficulty  in  evaluating  B0  and  by  analogy  with 
heat  transfer  a semiempirical,  semitheoretical  relationship  is  set  up 
for  mass  transfer.  Thus  it  is  customarily  considered  on  a semiempirical 
basis  that  the  mass  of  fuel  vapor  diffused  out  per  unit  time  under  the 
conditions  of  interest  here  is  given  by 


V = AoVfLa  U5) 

where  kg  is  the  mass  transfer  coefficient  for  the  film  in  the  case  of 
equimolal  diffusion.  Ranz  and  Marshall  (ref.  l)  present  a correlation 
for  kg,  but  again  it  was  necessary  to  extrapolate  their  correlation  to 
conditions  far  removed  from  the  experimental  conditions  under  which  it 
was  obtained.  In  making  this  extrapolation  the  line  of  reasoning  was 
as  follows:  The  correlation  was  for  the  film.  Therefore,  when  the 

density  of  the  air  and  of  the  film  were  markedly  different,  the  mean 
density  of  the  film  pm  should  be  used  instead  of  p&.  Another  choice 

to  be  made  in  the  extrapolation  was  whether  the  Ranz  and  Marshall  corre- 
lation is  for  kg  or  for  the  product  of  kg  and  a,  since  the  data  for 
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the  Ranz  and  Marshall  correlation  were  taken  under  conditions  where  a 
was  essentially  unity  and  therefore  not  a variable.  The  factor  pa  is 

sometimes  used  as  a correction  factor  instead  of  a (ref.  4).  However, 
when  it  is  so  used,  it  is  a logarithmic  mean  average,  while  in  all  of 
their  calculations  for  their  correlations  Ranz  and  Marshall  used  an 
arithmetic  mean  average.  The  difference  between  the  two  averages  is 
almost  undetectable  over  the  range  of  conditions  covered  in  their  experi- 
ments and  thus  it  is  not  clear  whether  their  correlation  included  a 
via  the  factor  pa.  It  is  also  not  clear  whether  the  use  of  the  arith- 
metic mean  rather  than  the  logarithmic  mean  is  an  essential  part  of  their 
correlation.  For  extrapolation  purposes  for  these  computations,  it  was 
decided  to  use  pa  as  an  arithmetic  average  and  to  assume  that  a was 

not  included  in  the  correlation  of  Ranz  and  Marshall.  Thus  the  final 
equation  used  was 

%u'  - - 2 + °-6(%o)l/3(HRe)l/2  <l6> 


The  factors  in  the  mass  transfer  equations  were  also  computed  for 
the  average  film  conditions  (see  appendix  B).  The  values  of  the  dif- 
fusion coefficient  Dy  were  arrived  at  by  using  the  techniques  recom- 
mended by  Hirschfelder,  Curtiss,  and  Bird  (ref.  5)*  It  can  be  seen  again 
from  the  above  equation  that,  to  calculate  the  mass  transfer  at  any 
instant,  a knowledge  of  the  heat  transfer  and  of  the  radius  and  velocity 
of  the  droplet  at  that  instant  is  essential. 


Velocity  Equation 

The  velocity  of  the  droplet  at  any  instant  during  the  vaporization 
process  is  essential  to  the  calculation  of  both  the  heat  and  mass  trans- 
fer at  that  instant,  since  it  enters  into  both  correlations.  The  droplet 
is  assumed  to  be  injected  with  a certain  initial  velocity  relative  to 
the  air.  Aerodynamic  drag  forces  will  then  either  slow  it  down  or  speed 
it  up  so  that  its  velocity  approaches  that  of  the  air.  The  less  the 
density  of  the  air,  that  is,  the  lower  its  pressure  and  the  higher  its 
temperature,  the  less  effective  the  drag  forces  and  the  greater  and  more 
rapid  the  penetration  through  the  combustion  chamber. 

From  aerodynamic  theory  the  drag  force  exerted  by  a moving  fluid 
upon  an  immersed  body  may  be  calculated  from  the  equation 


F = d>Acpm  X- 


(17) 
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where  Cd  is  the  coefficient  of  drag  for  the  given  body  form  and  is  a 
function  of  the  Reynolds  number;  Ac,  the  cross-sectional  area  of  the 
body;  and  pm,  the  density  of  the  fluid. 

For  a spherical  droplet: 


F = 


or 


dU  _ _ 3 cD 

dO  ~ 8 D Pl  rQ 


(18) 


Here  Cp  is  taken  as  the  drag  coefficient  for  spheres  (ref.  6)  at 
the  Reynolds  number  for  which  the  calculation  is  performed.  It  has  been 
pointed  out  by  Hughes  and  Gilliland  (ref.  7)  that  the  drag  coefficients 
for  solid  spheres  and  for  droplets  are  not  equal.  However,  the  condi- 
tions for  which  the  corrections  are  shown  do  not  correspond  to  the  con- 
ditions for  which  the  calculations  were  performed  and  therefore  no 
correction  was  made.  Furthermore,  extrapolation  of  the  data  of  Hughes 
and  Gilliland  would  indicate  that  the  correction  should  be  close  to 
unity  for  the  conditions  considered  here.  It  is  again  evident  that,  to 
calculate  the  deceleration  of  the  droplet  at  any  instant,  a knowledge 
of  both  the  heat  and  mass  transfers  is  essential. 


Method  of  Calculation 

The  dependency  of  the  heat  transfer,  mass  transfer,  and  velocity 
equations  upon  each  other  and  the  complexity  of  the  equations  neces- 
sitated a stepwise  integration  process  where  the  values  required  are 
evaluated  for  small  increments  of  droplet  temperature  over  a small 
period  of  time  during  which  all  properties  of  the  droplet  and  fluid  may 
be  assumed  to  be  constant.  A sample  calculation  is  shown  in  appendix  A. 
The  procedure  followed  was  first  to  assume  an  increment  of  time  £B,  a 
final  velocity  Ug,  and  a final  radius  r02  for  a certain  increment  of 

droplet  temperature  ^£Pl,  the  initial  values  of  velocity  TJq,  radius  rQq, 
and  temperature  Tq  for  that  increment  being  already  known  from  the 
previous  step,  or,  if  this  were  the  first  increment,  being  specified. 

For  that  increment  the  mean  values  of  the  droplet  temperature  T^, 
radius  r0,  and  velocity  U are  calculated.  Knowing  the  pressure  and 

temperature  of  the  air,  which  are  assumed  to  remain  constant  throughout 
the  process,  the  various  properties  in  the  film  are  calculated  (see 
appendix  B) . 
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Then  £U/£0  is  calculated  from  equation  (18)  and  U2  is  evaluated: 


«2  “ "1  - & (19) 

This  is  repeated  several  times,  correcting  for  Cp  each  time,  until 
a value  of  U2  is  found  that  satisfies  the  velocity  equation  and  the 
Reynolds  number  is  fixed  at  the  mean  velocity  of  the  increment  U. 

The  next  step  would  be  to  calculate  the  mass  transfer  coefficient  kg 

and  the  average  mass  transfer  w in  pounds  per  second  for  the  increment 
from  equations  (15)  and  (l6).  This  requires  a knowledge  of  the  partial 
pressure  of  the  fuel.  The  total  mass  transfer  for  the  increment  is  then 
found  from  w and  the  average  droplet  mass  m would  be 


m^_  + m2 

m 

2 


ml  + (ml  “ w 
2 


(20) 


In  this  equation  m^  would  be  specified  or  known  from  calculations  of 
the  previous  step. 

The  heat  transfer  to  the  droplet  surface  CL  is  then  calculated 

from  equations  (5)  and  (4)  and  the  remaining  calculations  take  place  as 
follows : 


- Qv  - w-^ 

For  all  calculations  shown  here  it  was  assumed  that  the  thermal  conduc- 
tivity of  the  fuel  was  high  and  the  droplet  temperature  was  uniform.  The 
validity  of  these  assumptions  is  discussed  in  appendix  C.  Therefore, 


A8 


ro2  = 


(21) 


(22) 


1 


\ 


NACA  TN  5179 


17 


The  values  of  Ug,  A0,  and  rQ2  calculated  from  equations  (19) , 

(21),  and  (22)  are  compared  with  those  assumed  at  the  beginning  of  the 
computations  and  are  usually  recalculated  at  least  once  more.  The  final 
values  for  this  increment  then  become  the  initial  values  of  the  next 
increment  and  the  computations  are  continued  for  small  increments  of 
ATl  until  Ql  is  found  to  approach  zero.  At  this  point  the  temperature 

asymptote  is  about  to  be  reached.  Here  the  computations  were  stopped, 
since  only  the  unsteady  state  was  being  studied.  Since  the  asymptotic 
temperature  was  not  known  for  most  conditions  prior  to  each  calcula- 
tion, ATl  and  the  number  of  steps  were  difficult  to  determine.  It  was 
found  that  between  12  and  16  steps  usually  gave  sufficient  accuracy. 

A model  1,  IBM  Card-Programmed  Electronic  Calculator  was  used  for 
the  computations.  The  special  program  to  solve  this  problem  uses  the 
floating-decimal  general-purpose  connections  of  the  calculator . One 
complete  calculation  for  specified  droplet  conditions  of  initial  radius, 
temperature,  and  velocity  and  air  conditions  of  pressure  and  temperature 
takes  approximately  2.0  hours  on  the  CPC  not  counting  the  setup  time, 
computation  of  the  properties  at  the  successive  droplet  temperature  incre- 
ments, and  card  prepunching. 

A typical  vaporization  history  calculated  by  this  method  is  shown 
in  figure  2 for  a droplet  of  n-octane  of  50-micron  initial  radius, 

50°  F initial  temperature,  and  100-fps  initial  velocity  with  respect  to 
air,  vaporizing  under  air  ambient  conditions  of  0.5-atmosphere  pressure 
and  1,000°  F temperature. 


Steady  State 

The  droplet  temperature  approaches  the  wet-bulb  temperature  asymp- 
totically during  the  vaporization  process.  At  the  wet -bulb  or  asymptotic 
temperature 


Qv  = Qa  (25) 

Substituting  values  from  equation  (5)  and  changing  TL  to  Tas  there 
is  obtained: 


Mo(%  - Tas)(^77)  - “X 


(24) 
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Substituting  vrCp^hA0  for  z and  rearranging: 

T = T- — (ez  - 1)  (25) 

as  B 0pf 

This  is  a simple  equation  for  the  asynqjtotic  or  wet-bulb  temperature 
as  a function  of  the  air  temperature  and  the  mass  and  heat  transfers.  A 
knowledge  of  the  latter  two  values  is  essential  for  a solution  of  this 
equation.  It  can,  however,  be  rewritten  using  only  the  properties  of 
the  liquid  and  film  as  follows:  Combining  equations  (2k)  and  (15)  and 

rearranging  give: 


T 

B 


T 


as 


PfLX 


a ^g 

z h 


e 


z 


1 


Combining  the  above  equation  with  equations  (5)  and  (l6)  and  rearranging: 


Tas  “ tB 


^yPmPfL* 
Pa^m 


a 


Ni 


Nu 


N, 


(26) 


Nu 


ez  - 1 


Solution  of  equation  (26)  must  necessarily  be  by  trial  and  error, 
since  the  component  properties  are  all  functions  of  both  Tag  and  Tg 

as  well  as  of  the  air  pressure  and  the  fuel  used.  The  factor  z/(ez  - l) 
must  be  approximated  since  this  factor  is  affected  by  the  volatility  of 
the  fuel.  It  can  be  stated  that  for  very  low  volatility  fuels  this  factor 
can  usually  be  taken  as  unity  with  little  error.  For  high- volatility 
fuels  at  low  air  temperatures  and  high  air  pressures  z/(ez  - l)  can 
still  be  taken  as  unity  with  little  error.  For  example,  using  n-pentane, 
this  factor  has  a value  of  0.957  at  100°  F air  temperature  and  5-atmosphere 
air  pressure  while  n-octane,  a less  volatile  fuel,  has  a value  of  0.998 
for  this  factor  at  the  same  conditions.  However,  for  high-volatility 
fuels  at  high  air  temperatures  and  low  air  pressures  lower  values  have  to 
be  used.  n-Pentane  has  a value  of  about  O.65  at  500°  F and  1 atmosphere 
air  conditions. 


The  mass  transfer  during  the  steady-state  period  of  vaporization  can 
be  determined  by  combining  equations  (15)  and  (16)  as  follows: 


2DvPmPfLa  C 
rn  2 + 

D O 


W 
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Since  the  temperatures  across  the  film  have  become  essentially  constant, 
the  above  equation  can  be  rewritten  as: 


» - V.|f  + 


(27) 


where  Kp  and  Kq  are  constants. 


Equation  (27)  shows  that  for  a particular  fuel  the  rate  of  mass 
transfer  during  the  steady  state  is  a function  of  only  the  radius  and 
velocity  of  the  droplet  at  any  instant,  provided  the  air  temperature 
and  pressure  remain  unchanged. 


Equation  (26)  for  the  steady  state  where  the  temperatures  across 
the  film  are  fixed,  and  therefore  all  the  properties  in  it,  can  be 
rewritten  as 


T 


as 


= T- 


■B 


K 2 - Mrou)1/2 
3 2 + K^u)1/2 


(28) 


where  Kj,  K^,  and  are  constants  for  any  particular  fuel  and  air 

conditions  at  the  asymptotic  or  wet-bulb  temperature  Tas.  Therefore, 

Tas  depends  upon  the  ratio  of  the  Nusselt  numbers  for  mass  and  heat 

transfer.  It  can  be  deduced,  after  examining  the  above  equation,  that 
the  wet-bulb  temperature  is  dependent  on  the  droplet  radius  and  velocity. 
At  high  droplet  speeds  or  for  relatively  large  droplet  radii  the  fac- 
tor 2 in  both  the  Nusselt  numbers  becomes  relatively  negligible  and  the 
ratio  of  the  Nusselt  numbers  and  consequently  the  wet-bulb  temperature 
become  essentially  constant.  However,  at  lower  values  of  droplet  speed 
or  droplet  radius  the  wet-bulb  temperature  rises  slightly  until  either 
the  speeds  or  radii  approach  zero,  where  the  Nusselt  number  ratio 
approaches  unity  and  the  wet-bulb  temperature  is  again  fixed. 

This  dependence  proved  to  be  small  for  most  conditions  of  interest 
in  this  paper.  Figure  3(a)  shows  ATas  as  calculated  for  various  veloci- 
ties, where  zXEas  is  defined  as  the  deviation  of  the  actual  wet-bulb  tem- 
perature from  that  calculated  at  large  Reynolds  numbers*  figure  3(b) 
shows  ATas  for  various  radii  for  n-pentane  at  1-atmosphere  air  pressure 
and  100°  F air  temperature.  It  should  be  noted  that  no  attempt  was  made 
to  include  a dependence  of  the  vapor  pressure  of  fuels  on  droplet  size 
as  the  droplet  radius  became  very  small  in  calculating  the  information 
shown  in  figure  3(t>) • It  is  to  be  expected  that  the  rise  of  vapor  pres- 
sure as  corrected  for  very  small  droplet  radii  would  lower  the  values  of 
the  wet-bulb  temperature  again  at  small  drop  radii. 


20 


NACA  TN  3179 


Further  examination  of  equation  (26)  reveals  that  the  wet-bulb  tem- 
perature approaches  the  air  temperature  in  case  of  low-volatility  fuels. 
This  is  to  be  expected,  since  the  mass  transfer  has  to  be  high  enough  to 
carry  away  the  entire  heat  transferred  to  the  droplet.  In  the  case  of 
no  mass  transfer,  the  wet-bulb  temperature  equals  the  air  temperature. 


APPARATUS 


After  investigating  theoretically  the  importance  of  the  unsteady 
state  during  the  vaporization  period  and  conducting  a number  of  calcu- 
lations in  the  unsteady  state,  it  was  felt  that  experimental  verifica- 
tion of  its  importance  was  desirable.  Accordingly,  experimental  appara- 
tus has  been  set  up  in  an  attempt  to  establish  the  importance  of  the 
unsteady  state  and  a few  experimental  records  have  been  obtained. 


Figure  4 shows  a schematic  diagram  of  the  apparatus  constructed  to 
determine  the  temperature -time  histories  of  vaporizing  droplets.  Air 
is  tapped  from  the  house  line  at  100  Ib/sq  in.  abs  and  passed  through 
a porous-stone  filter  to  remove  water,  oil,  and  other  impurities.  Its 
pressure  is  then  controlled  by  a pressure  regulator,  after  which  it 


passes  through  a 2—  -inch  pipe  where  its  pressure  and  temperature  are 


measured  by  a precision  Bourdon  gage  and  a thermometer.  The  diameter  of 
the  pipe  at  this  measuring  section  is  made  large  enough  so  that  the  dif- 
ference between  total  and  static  temperatures  is  negligible  at  the  highest 
rate  of  flow  to  be  used.  The  air  is  then  passed  through  either  the 
3/32-  or  the  3/16- inch-diameter  critical-flow  orifice.  The  air  is  then 
heated  with  an  electrical  resistance  heater,  after  which  it  passes  through 
a calming  section  to  a 14- inch- long  section  of  2-inch  flexible  tubing 
which  terminates  in  a 1-inch  I.S.A.  standard  flow  nozzle.  The  exit  sec- 
tion of  the  nozzle  is  covered  with  a l60-mesh  screen.  The  screen  assists 
in  equalizing  the  velocity  over  the  nozzle  opening  to  a value  calculated 
from  the  average  volumetric  velocity. 


The  liquid  droplet  is  hung  on  the  junction  of  the  thermocouple 
about  l/2  inch  above  the  exit  section  of  the  nozzle.  A 0.010- inch  iron- 
constantan  thermocouple  was  used  for  the  preliminary  investigations 
reported  herein.  A clean  and  small  thermocouple  junction  was  found  essen- 
tial for  a minimum  of  droplet  distortion.  It  was  found  that  droplets  of 
700-  to  900-micron  radius  were  the  best  size  to  be  used  with  the  above 
thermocouple  wire  size  without  distortion  during  the  major  portion  of  the 
heating-up  and  vaporization  period. 

An  illuminating  system  was  provided  for  viewing  the  droplet.  A 
300-watt  General  Electric  T8|  projection  lamp  was  used.  A heat-absorbing 
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glass  was  placed  in  the  path  of  the  light  to  minimize  the  heating  of 
the  droplet  by  radiation  from  the  lamp.  An  optical  system  having  a maxi- 
mum magnification  of  86  was  provided  for  measuring  the  diameter  of  the 
droplet.  The  optical  system  was  calibrated  by  placing  a wire  of  known 
diameter  in  place  of  the  droplet. 

The  experimental  procedure  was  as  follows:  The  air  was  heated  to 

the  desired  temperature  with  the  exit  nozzle  swung  away  from  the 
thermocouple . 

The  exit  velocity  of  the  air  stream  was  adjusted  by  varying  the 
inlet  air  pressure  ahead  of  the  critical  flow  orifice.  Low  air  veloci- 
ties were  used  to  avoid  blowing  the  droplet  off  the  thermocouple  junction. 

The  droplet  was  then  formed  and  suspended  on  the  thermocouple  using 
a syringe  and  hypodermic  needle.  The  initial  droplet  radius  was  then 
measured  by  the  optical  magnification  system.  The  air  nozzle  at  the 
end  of  the  flexible  tubing  section  was  then  suddenly  swung  under  the 
droplet.  The  output  of  the  thermocouple  was  measured  and  recorded  by 
a 1-second  full-scale  model  G Leeds  & Northrup  electronic  self-balancing 
Speedomax  potentiometer.  The  temperature  of  the  droplet  was  measured  by 
the  calibrated  thermocouple  output,  while  the  time  scale  was  computed 
from  the  known  speed  of  the  recorder  chart. 


RESULTS  AND  DISCUSSION 


The  typical  calculated  vaporization-time  histories  shown  in  fig- 
ures 2 and  5 through  8 contain  the  following  information: 

(1)  Droplet  temperature,  T^,  °F 

(2)  Droplet  mass,  percent  of  original,  M,  percent 

(3)  Radius,  rQ,  inches 

(4)  Velocity  relative  to  the  air,  U,  fps 

(5)  Penetration  or  distance  traveled  relative  to  the  air,  d,  inches 

(6)  Instantaneous  rate  of  mass  transfer,  w,  pounds  per  second 

(7)  Instantaneous  total  heat  transfer  from  the  air  to  the  film,  Q, 

Btu/sec 

(8)  Instantaneous  heat  of  vaporization,  Q^,  Btu/sec 

(9)  Instantaneous  heat  transfer  to  the  liquid  droplet,  Q^,  Btu/sec 

The  fuel,  initial  droplet  radius,  initial  droplet  temperature,  initial 
droplet  speed  relative  to  the  air,  air  pressure,  and  air  temperature  are 
indicated  in  the  title  of  each  figure.  The  portion  of  the  total  heat  trans- 
fer Q that  is  carried  back  with  the  diffusing  vapor  in  the  form  of 
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superheat  can  he  determined,  from  the  curves  by  taking  the  difference 
between  Q and  the  sum  of  and  Q^. 

Figures  2 and  5 through  8 represent  the  results  of  only  7 out  of 
the  252  cases  for  which  calculations  were  performed.  A complete  set  of 
graphs  presenting  the  results  of  all  the  calculations  is  contained  in 
the  manuscript  copy  of  this  report.  This  copy  is  available  for  loan  or 
reference  in  the  Division  of  Research  Information,  National  Advisory 
Committee  for  Aeronautics,  Washington,  D.  C. 

Examination  of  these  typical  graphs  (figs.  2 and  5 to  8)  yields  the 
following  conclusions: 

For  low- volatility  fuels,  for  example,  figure  2 for  n-octane, 
speaking  of  the  heat  transferred  to  the  surface  of  tne  liquid  droplet, 
the  amount  of  heat  going  to  sensible  heat  far  exceeds  the  amount  of  heat 
going  to  supply  the  latent  heat  of  vaporization  during  most  of  the  unsteady- 
state  or  heating-up  period  and  especially  right  after  injection.  The 
higher  the  air  temperature  or  the  higher  the  air  pressure,  the  greater  the 
portion  going  to  sensible  heat.  Only  as  the  droplet  approaches  the  steady 
state  does  the  relation  reverse  itself,  as  required  by  the  fact  that  Ql 

becomes  negligible  and  essentially  equal  to  zero  during  the  steady  state. 
This  effect  is,  of  course,  largely  dependent  on  the  initial  injection 
temperature  of  the  liquid  fuel.  At  high  initial  liquid  temperatures  sen- 
sible heat  can  be  lost  by  the  droplet  or  a net  heat  transfer  Q from  the 
liquid  to  the  air  may  even  occur  (e.g.,  fig.  5 for  n-octane  at  200°  F 
initial  liquid  temperature  and  100°  F air  temperature) . 

For  high- volatility  fuels,  however  (e.g.,  fig.  6(a)  for  n-pentane 
at  1-atmosphere  air  pressure  and  500°  F air  temperature),  the  latent 
heat  of  vaporization  exceeds  the  amount  going  to  sensible  heat.  In  some 
cases,  especially  at  low  air  pressures  (e.g.,  fig.  6(b)  for  n-pentane), 
the  liquid  undergoes  a cooling  down  or  loses  sensible  heat  at  0.5-atmosphere 
air  pressure  regardless  of  the  air  temperature.  In  this  case  the  net  heat 
transfer  from  the  air  goes  entirely  into  latent  heat  of  vaporization.  The 
balance  of  the  heat  required  to  supply  the  latent  heat  of  vaporization  is 
drawn  from  the  droplet  itself.  Only  at  high  air  temperatures  and  high  air 
pressures  does  the  sensible -heat  requirement  of  the  droplet  exceed  that  of 
the  heat  of  vaporization,  and  only  for  a shorter  portion  of  the  unsteady 
state . 

The  rate  of  mass  transfer  usually  increases  with  time  as  the  droplet 
is  heating  up.  The  increase  is  much  more  pronounced  at  higher  air  tem- 
peratures and  lower  air  pressures  (e.g.,  fig.  2 for  n-octane  at 
0.5-atmosphere  air  pressure  and  1,000°  F air  temperature).  At  high  air 
pressures  and  low  air  temperatures  (e.g.,  figs.  7 and  8(a)  for  benzene 
and  n-octane  at  5-atmosphere  air  pressure  and  100°  F air  temperature) 
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the  rate  of  mass  transfer  undergoes  a decrease,  although  the  droplet  is 
heating  up.  The  reason  for  this  is  that  in  this  case  the  denser  air 
slows  down  the  droplet  very  rapidly  with  a consequent  rapid  drop  in  the 
Reynolds  number  which  offsets  the  favorable  effect  of  the  exponential 
increase  in  the  vapor  pressure  of  the  liquid  fuel  with  temperature  on 
the  rate  of  mass  transfer  w.  At  intermediate  cases,  a rise  in  w 
followed  by  a drop  can  be  noticed  (e.g.,  fig.  8(b)  for  n-octane  at 
1-atmosphere  air  pressure  and  100°  F air  temperature).  In  the  case  of 
droplets  cooling  down  during  the  unsteady  state,  the  rate  of  mass  trans- 
fer consistently  decreases  with  time.  In  the  steady  state  the  rate  of 
mass  transfer  becomes  a function  of  only  the  radius  and  velocity  of  the 
droplet  and  thus  undergoes  a small  but  steady  decrease  in  value  (see 
eq.  (27)). 

For  comparison  purposes  the  three-dimensional  graphs  shown  in 
figs.  9 through  12  are  presented.  These  are  temperature -time  history 
curves  for  the  three  fuels  at  different  fuel  and  ambient  air  conditions. 

In  all  these  curves  the  initial  temperature  of  the  droplets  was  5O0  F 
and  the  initial  velocity  relative  to  the  air  was  100  fps. 

Figure  9(a)  shows  temperature -time  histories  of  n-octane  at  an  air 
temperature  of  500°  F,  for  various  droplet  sizes  and  air  pressures.  It 
can  be  seen  that  for  the  same  air  pressures  but  different  radius  droplets 
the  asymptotic  temperatures  are  equal  and  the  three  curves  representing 
the  three  sizes  at  any  one  air  pressure  form  a horizontal  plane  after  all 
three  droplets  have  reached  their  steady  states.  The  ends  of  the  hori- 
zontal planes  shown  by  the  broken  lines  have  no  particular  meaning  and 
were  drawn  for  illustrative  purposes  only,  since  the  10-micron  droplets 
would  vaporize  before  reaching  this  point  and  the  100-micron  droplets 
would  last  longer.  Figure  9(b)  shows  similar  temperature-time  histories 
for  n-octane  except  that  the  air  pressure  was  kept  constant  at  1 atmos- 
phere and  the  air  temperature  varied.  Figure  9(c)  shows  temperature- 
time histories  for  n-octane  where  all  the  droplets  had  a 50-micron 
initial  radius  and  both  the  air  temperature  and  pressure  were  varied. 

In  this  case  none  of  the  nine  temperature -time  histories  shown  would 
have  the  same  asymptotic  temperature.  The  end  lines  joining  the  curves 
at  the  same  air  temperature  represent,  therefore,  a variation  of  asymp- 
totic temperatures  with  air  pressures  at  three  different  air  temperatures. 

Figures  10(a)  to  10(c)  and  ll(a)  to  ll(c)  show  the  same  histories 
described  in  figures  9(a)  to  9(c)  but  for  benzene  and  n-pentane, 
respectively. 

Figure  12(a)  shows  temperature -time  histories  of  the  three  fuels 
where  the  initial  droplet  radius  was  50  microns,  the  air  pressure  was 
1 atmosphere,  and  the  air  temperature  was  varied.  Figure  12(b)  shows 
the  three  fuel  droplets  having  the  same  initial  radius,  with  the  air 
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temperature  kept  constant  at  500°  F and  the  air  pressure  varied.  It 
can  be  seen  from  these  curves  that  the  lower  volatility  fuels  have 
higher  asymptotic  temperatures  and  take  a longer  time  to  reach  them  for 
the  same  operating  and  ambient  conditions.  By  plotting  the  results  of 
the  present  calculations,  it  was  found  that,  with  all  other  conditions 
constant,  the  time  for  droplets  to  reach  the  wet-bulb  temperature  varied 
approximately  as  the  initial  radius  to  the  1.8  power. 

In  the  combustion  chambers  the  mass  of  the  fuel  droplets  vaporized 
and  the  penetrations  of  these  droplets  relative  to  the  air  are  essential 
in  deciding  the  zone  where  combustion  is  likely  to  take  place.  Fig- 
ures 13(a)  through  13(i)  show  time  histories  of  droplet  penetration  d, 
relative  to  the  air,  and  mass,  percent  of  original,  M during  the 
unsteady  state  for  n-octane  at  nine  different  ambient  air  conditions. 
Each  of  these  figures  contains  the  previous  information  for  the  three 
initial  droplet  sizes  and  two  initial  droplet  velocities.  As  the  tem- 
perature decreases  but  with  the  same  air  pressure,  the  droplets  take 
a longer  time  and  penetrate  farther  before  reaching  their  steady  states j 
as  the  air  pressure  increases  but  for  the  same  air  temperature,  the 
droplets  take  a longer  time  but  penetrate  less  before  reaching  their 
steady  states. 

Figure  lh-  presents  a comparison  of  benzene  and  n-octane  wet-bulb 
temperatures  calculated  in  this  report  with  wet-bulb  temperatures 
experimentally  obtained  by  Ingebo  (ref.  2).  Referring  to  the  benzene 
curves , it  can  be  seen  that  the  calculated  wet-bulb  temperatures  are 
always  somewhat  higher  than  the  experimental  wet-bulb  temperatures. 

This  would  indicate  that  the  computed  mass  transfer  was  too  low  or  the 
computed  heat  transfer  was  too  high  or  both.  Referring  to  the  n-octane 
curves,  it  can  be  seen  that  there  is  a cross-over  of  the  two  curves. 
Extrapolation  of  the  data  of  Ingebo  would  indicate  that  the  droplet 
would  boil,  in  contradiction  to  the  diffusion  theory.  Whether  this  is 
a deficiency  in  the  theory  or  in  the  experimental  data  of  Ingebo  is  not 
clear. 


A typical  temperature -time  history  obtained  in  the  experimental 
investigations  is  shown  in  figure  15  for  a droplet  of  n-octane  of  approxi- 
mately 750-micron  initial  radius  when  exposed  to  air  moving  at  a constant 
velocity  of  about  l6  fps.  The  portion  1-2  of  this  graph  shows  the  output 
of  the  thermocouple  when  initially  in  room  air.  At  point  2 the  droplet 
was  hung  on  the  thermocouple.  Portion  2-3-4  of  the  graph  represents  the 
rapid  response  of  the  thermocouple  to  the  cooler-than-room-air  droplet 
followed  by  a gradual  cooling  to  the  wet-bulb  temperature  corresponding 
to  room  air  at  point  4.  This  period  is  of  particular  interest  in  that 
it  gives  an  indication  of  the  speed  of  response  of  the  thermocouple  and 
recording  instrument. 

At  point  b a stream  of  hot  air  was  suddenly  swung  under  the  droplet. 
The  unsteady-state  or  heating-up  portion  of  the  curve  can  be  seen  from 
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point  4 to  a point  chosen  arbitrarily  at  5 where  the  droplet  would  be 
Judged  to  have  closely  reached  the  wet-bulb  temperature  corresponding  to 
the  temperature  of  the  air  stream.  The  portion  5-6  of  the  graph  is  con- 
sidered to  be  the  steady-state  vaporization  period.  The  last  part  of 
this  portion,  near  6,  is  of  questionable  reliability,  since  observation 
showed  that  as  the  droplet  nears  complete  vaporization  it  deviates  from 
the  spherical  shape  assumed  in  the  calculations  and,  instead,  takes  a 
distorted  shape  depending  on  the  exact  shape  of  the  thermocouple  Junc- 
tion. It  should  also  be  noted  that  the  purity  of  the  fuels  was  only 
about  98  percent  and  that  this  portion  of  the  curve  may  represent  the 
vaporization  of  the  high-boiling -point  impurities  in  the  mixture. 

At  point  6 a thin  film  of  the  liquid  fuel  still  remains  on  the 
thermocouple  element.  This  soon  disappears  and  the  thermocouple  tem- 
perature rises  to  the  air  temperature  as  represented  by  point  7*  At 
that  point  the  thermocouple  is  completely  dry  and  the  temperature 
recorded  at  7 is  the  air  temperature  during  the  entire  plot.  Since 
this  takes  only  a short  time,  the  air  temperature  is  not  likely  to  have 
varied  during  this  period.  This  procedure  has  the  advantage  of  using 
a single  thermocouple  to  measure  both  the  droplet  and  air  temperatures . 

Since  the  amount  of  liquid  fuel  remaining  as  a thin  film  on  the 
thermocouple  element  at  point  6 in  the  graph  is  minute  compared  with 
the  total  liquid  contained  in  the  initial  droplet,  the  relative  magni- 
tudes of  time  of  portions  4-5  and.  5-8  of  the  graph  can  be  assumed  to 
represent  the  relative  importance  of  the  unsteady  and  steady  states, 
respectively. 

Figure  16  shows  several  time  histories  of  droplets  of  n-octane  all 
of  approximately  equal  size  as  subjected  to  air  of  varying  temperature. 
The  portions  of  the  curve  corresponding  to  2-3-4  in  figure  15  have  been 
omitted  here.  It  can  be  seen  from  these  curves  that  the  unsteady-state 
portion  of  the  vaporization  time  assumes  greater  importance  and  becomes 
dominant  at  higher  air  temperatures. 

Figure  17  shows  time  histories  of  liquids  of  varying  volatility 
vaporizing  under  the  same  conditions  of  pressure,  temperature,  and 
velocity. 

Since  in  the  calculated  curves  presented  in  this  report  the  velocity 
of  the  droplet  was  made  to  decrease  relative  to  the  air  because  of  aero- 
dynamic drag  forces,  while  the  experimental  curves  were  measured  at 
essentially  constant  air  velocity  past  the  droplets,  caution  must  be 
exercised  in  comparing  the  experimental  results  thus  obtained  with  the 
calculated  curves.  Figure  18  shows  two  calculated  temperature,  velocity, 
and  rate  of  mass  transfer  (unsteady -state)  histories  of  two  droplets  of 
n-octane  of  50-micron  initial  radius  and  100-fps  initial  velocity  each. 
One  was  made  to  slow  down,  while  the  other  remained  at  100-fps  velocity 
throughout  the  calculated  period.  Other  conditions  were  kept  equal. 
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It  can  be  seen  that  the  constant-velocity  droplet  spends  less  time  in 
the  unsteady-state  period  but  has  a higher  rate  of  mass  transfer  during 
that  period  with  consequently  shorter  total  vaporization  time.  The  net 
result  would  be  that  the  unsteady  state  is  relatively  longer  for  the 
droplet  that  is  moving  at  a constant  speed. 

The  following  deductions  can  now  be  made. 

(1)  For  the  same  fuel  the  unsteady  state  becomes  relatively  larger 
in  magnitude  with  respect  to  the  total  vaporization  time  the  higher  the 
air  temperature. 

(2)  For  the  same  air  temperatures  the  unsteady  state  is  relatively 
larger  in  magnitude  with  respect  to  the  total  vaporization  time  the  higher 
the  volatility  of  the  fuel. 

(3)  High- volatility  fuels  have  lower  wet-bulb  temperatures  for  the 
same  air  conditions. 

(4)  At  relatively  low  air  temperatures  extremely  low  volatility 
fuels  have  wet -bulb  temperatures  close  to  the  air  temperatures  and  spend 
only  a very  small  portion  of  their  vaporization  time  in  the  unsteady 
state.  At  high  air  temperatures,  however,  the  difference  between  the 
air  and  wet-bulb  temperatures  increases. 

(5)  For  any  one  fuel  the  wet-bulb  temperature  is  higher  the  higher 
the  air  pressure  or  temperature. 

The  penetration  - mass -vaporized  relationships  of  different  size 
droplets  are  of  particular  interest  in  jet-engine  combustion-chamber 
design.  When  the  fuel  is  injected  continuously  into  a stream  of  moving 
air,  it  leaves  the  injector  nozzle  first  as  a ligament  or  sheet  after 
which  this  ligament  or  sheet  breaks  down  to  different  size  droplets 
originally  moving  at  the  same  speed . If  the  velocity  of  these  droplets 
is  higher  than  the  air  velocity,  these  droplets  start  slowing  down  rela- 
tive to  the  air,  with  the  smaller  droplets  slowing  down  faster,  and 
penetrate  a lesser  distance  while  losing  their  mass  faster  and  earlier 
along  their  travel  path.  The  mass  of  vapor  given  away  by  each  droplet 
is  carried  along  with  the  air  through  the  combustion  chamber.  To  this 
vapor  is  continuously  added  the  vapor  given  away  by  the  larger  droplets 
that  have  been  formed  later  in  time  but  that  reached  that  point  because 
of  their  more  rapid  and  greater  penetration.  The  total  amount  of  fuel 
vapor  present  at  a certain  cross  section  of  the  combustion  chamber  there- 
fore consists  of  the  total  vapor  given  away  by  all  the  drops  before 
reaching  that  cross  section.  In  order  to  have  a short  combustion  cham- 
ber, a combustible  air-vapor  mixture  should  be  formed  a short  distance 
after  injection.  This  would  necessitate  that  the  largest  number  of  the 
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drop  sizes  fall  in  the  smaller  radius  range  so  that  the  weight  distri- 
bution of  different  droplet  sizes  in  the  spray  would  not  he  weighted  too 
highly  in  favor  of  the  droplets  of  large  radius  range,  since  these  large 
droplets  have  a much  greater  penetration  before  they  contribute  much 
vapor  toward  forming  a combustible  mixture. 

To  illustrate  the  possible  vise  of  the  curves  of  droplet  mass  versus 
time  for  combustion-chamber  design  purposes,  and  to  emphasize  the  impor- 
tance of  the  unsteady  portion  of  the  time  history  of  the  droplets,  the 
following  example  was  worked  out  with  certain  arbitrary  but  reasonable 
assumptions  made  as  to  the  couibustion-chamber  operating  conditions. 

These  assumptions  were: 

(1)  No  effect  of  one  evaporating  droplet  upon  the  other.  This  is 
in  obvious  contradiction  to  the  facts,  but  since  this  effect  is  unknown 
it  was  necessary  to  neglect  it. 

(2)  A constant  air  velocity  of  50  fps. 

(3)  An  air  pressure  of  1 atmosphere  and  a temperature  of  500°  F. 

The  500°  F air  temperature  is  admittedly  high  with  respect  to  the  air 
temperature  coming  from  the  compressor,  but  calculations  were  available 
for  this  temperature  and  recirculation,  which  was  neglected,  would 
undoubtedly  increase  the  air  temperature. 

(4)  The  fuel  used  was  n-octane. 

(5)  A drop-size  weight  distribution  of  25  percent  for  the  10-micron 
group,  65  percent  for  the  50-micron  group,  and  10  percent  for  the 
100-micron  group. 

(6)  An  initial  fuel  velocity  relative  to  the  air  of  100  fps. 

The  previously  performed  calculations  were  extended  through  the 
steady  state  until  the  droplets  nearly  disappeared/  Such  a calculation 
showing  temperature,  mass,  and  penetration  relative  to  the  air  versus 
time  for  the  50-micron  droplet  is  shown  in  figure  19.  The  masses  vapor- 
ized from  each  size  droplet  were  then  calculated  and  plotted  versus  dis- 
tance traveled  by  each  droplet  relative  to  the  combustion  chamber  using 
the  assumed  airspeed  of  50  fps,  as  shown  in  figure  20. 

Using  the  arbitrarily  chosen  droplet  size  distribution,  the  percent- 
age of  the  total  fuel  vaporized  from  each  droplet  size  group  up  to  any 
distance  after  droplet  formation  was  then  plotted  in  figure  21.  The 
total  vapor  present,  which  is  the  sum  of  the  three  curves,  is  also  shown 
in  figure  21.  Assuming  a fuel-air  mixture  strength  in  the  primary  com- 
bustion zone  to  be  in  the  neighborhood  of  0.08,  and  a possible  inflamma- 
bility ratio  of  0.03  of  the  vapor-air  ratio,  a total  vaporization  in  the 
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neighborhood  of  35  percent  would  be  adequate  for  initiating  combustion 
under  these  assumptions.  According  to  the  calculations  presented  in 
figure  21,  this  would  take  place  8 inches  after  the  droplets  have  been 
formed. 

It  is  to  be  noted  that  these  calculations  are  approximate  at  best. 
Turbulence,  backflow,  and  heat  transfer  by  radiation  have  not  been  taken 
into  account.  However,  they  serve  to  show  the  possible  use  of  the  cal- 
culated curves  for  estimation  purposes.  A knowledge  of  droplet  size 
distribution  seems  also  to  be  fundamental  in  any  such  vaporization 
calculations . 

An  important  observation  shown  by  the  above  calculations  is  that 
the  50-  and  100-micron  droplets  reach  the  combustion  zone  while  still 
in  the  heating-up  or  unsteady-state  portion  of  vaporization,  a fact 
that  further  emphasizes  the'  importance  of  the  unsteady  state  in  any 
combustion  studies.  It  is  the  opinion  of  the  authors  that  the  unsteady 
state  is  of  primary  importance  in  forming  a combustible  mixture  and  that 
experimental  studies  which  could  be  used  either  to  verify  or  to  correct 
the  extrapolation  of  the  heat  and  mass  transfer  correlations  used  should 
be  conducted.  Until  these  extrapolations  are  confirmed,  the  absolute 
accuracy  of  the  curves  can  be  questioned,  although  it  is  not  believed 
that  any  of  the  conclusions  reached  would  be  altered.  After  these  rela- 
tionships have  been  confirmed  or  corrected,  it  would  then  seem  necessary 
to  determine  the  effect  of  one  droplet  upon  another.  This  information 
together  with  drop  size  distribution  should  then  permit  a rational  and 
reasonable  determination  of  the  time  required  to  form  a combustible  mix- 
ture in  a heterogeneous  mixture  of  fuel  droplets  and  air. 


University  of  Wisconsin, 
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APPENDIX  A 

SAMPLE  OF  CALCULATIONS 


The  unsteady-state  portion  of  the  vaporization  time  of  a droplet  of 
benzene  was  calculated  for  the  following  data: 

Droplet  initial  radius  r •[_  = 50  microns  = I.9685  X 10“5  in. 

Droplet  initial  temperature  Ty_  = 50°  F = 510°  R 

Droplet  initial  velocity  relative  to  the  air  U]_  = 4 00  fps  = 

4,800  in. /sec 

Air  pressure  Pj  a 1 atm  = 14.696  lb/sq  in.  abs 

Air  temperature  Tg  = 1,000°  F = 1,460°  R 

A stepwise  method  was  used  in  integrating  the  equations  as  follows: 

The  temperature  of  the  droplet  was  assumed  to  undergo  a small  increase  ATl 
over  a small  period  of  time  A8,  during  which  the  properties  of  the  liquid 
fuel,  fuel  vapor,  and  air  could  be  considered  to  be  constant.  The  heat 
and  mass  transfer  during  this  period  were  then  calculated  together  with 
the  velocity  change,  from  which  the  time  required  for  the  temperature 
increase  together  with  the  radius  and  velocity  at  the  end  of  the  increment 
were  calculated.  However,  these  three  values  had  to  be  assumed  prior  to 
the  calculation  and  compared  with  the  calculated  values. 

For  a droplet-temperature  increase  of  5°  F it  was  assumed  that: 


A8  = 2.7  X 10-5  sec 


^2  - k>725  in. /sec 
ro2  = 1,970  x 10"5  in. 


For  this  increment: 


T-^  as  average  temperature  of  droplet  = 


510  + 515  o 

= 512.5  R 


2 
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The  various  properties  in  the  film  and  for  the  droplet  were  then 
calculated  and  read  as  follows: 

PfL  = 0.970  lb/sq.  in.  abs 
= 2.678  x 10"5  lb/cu  in. 

Pm  = 1.475  X 10"6  lb/(in.) (sec) 
kjjj  = 5.591  x 10-7  Btu/ ( in. ) ( sec)  ( °F) 

Cp^  = 0.2609  Btu/(lb)(°P) 

Dy  = 3.60  X 10"2  sq.  in. /sec 
A = 189.0  Btu/lb 
Cpf  = 0.422  Btu/(lb)(°F) 
as  3.154  x 10-2  lb/cu  in. 

0^2  = 3.150  x 10“2  lb/cu  in. 

= 0.4080  Btu/(lb)(°F) 

r0  = Average  radius  of  the  droplet  for  the  increment 

- rol  + ro2 
2 

_ ^1.9685 

1.9693  x 10-5  in. 
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U = Average  velocity  of  the  droplet 

U]i  + U2 

* 2 

- ^j800  + ^725 
2 

= 4,762.5  in. /sec 


%e  = fieynol(is  number 

2^Pm 

Mm 

= 340.568 


CD  (at  NRe)  = 0.6633 


Applying  equation  (l8): 


AO 

AB 


2 CD  — = 2.4326  in. /sec2 
8 PL1  ro 


and 


U = 65.681  in. /sec 


Therefore 


U2  = 4,800  - 65.681  = 4,734.319  in./sec 


as  compared  with  the  assumed  value  of  4,725  in./sec.  This  process  was 
then  repeated,  selecting  the  correct  coefficient  of  drag  each  time  until 
a value  for  Ug  was  arrived  at  that  satisfied  the  velocity  equation. 
Therefore 


U2  = 4,730.79  in./sec 
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%e  = 3^0.775 


and. 


Therefore 


NSc  = = 1.5299 


PmPv 


kg  ' S^lf  + °-6(%o)1/3(%e)l/^| 


= 2.5427  x lO-4  sec-1 


a = 


_T_ 

pfL 


= 1.030 


logg 


rT 


PT  “ pfL 


Therefore 


w = kgaPfL4rtr02  = 1.2381  x 10"8  lb/sec 


Cp  um 

Npr  = —2—  = 0.6883 


Therefore 


h = 


SjE  + °-6(SFr)1/5(»Re)l/3 


= I.6721  x 10“5  Btu/(sec)(sq.  in.)(^F) 


Q = hAQ(TB  - Tl)  = 7.7213  x 10-5  Btu/sec 


z 


e 


z 


1 


O.9683 
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Therefore 


and 


Qy  = Q - — = 7.4764  x 10" 5 Btu/sec 

ez  - 1 

Qs  = Q - Qv  = 2.4488  x 10"6  Btu/sec 
= wA  = 2.34006  x 10"6  Btu/sec 
Ql  = Qy  “ = 7*2424  x 10~5  Btu/sec 


= Mass  of  droplet  at  beginning  of  increment 
= 4/3«rol3PL  (at  510°  R) 

= 1.00936  X 10-9  lb 


m2  = Mass  of  droplet  at  end  of  increment 
= m1  - w £B 

= 1.00901  x 10-9  lb 

mi  + mo  Q 

m = — = = 1.009186  x 10"y  lb 

2 

Therefore 

mCpL 

A8  = = 2.8426  x 10"  3 sec 

<*L 


(as  compared  with  the  estimated  value  of  2.7  x 10"  ^ second). 
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3^ 


/ >2  V^5  v 

-02  - (j^j  - 1-W17  X 10-3  1„. 

(as  compared  with  the  estimated  value  of  1.970  x 10"5  inch). 

The  calculations  for  this  first  increment  were  then  repeated  taking 
the  calculated  values  of  time,  velocity,  and  radius  as  new  assumptions 
to  obtain  more  accurate  results.  The  final  values  of  U2,  r g,  an(*  m2> 

calculated  for  the  first  increment,  then  became  the  original  values  of 
the  next  increment  and  the  same  stepwise  method  was  used  for  the  second 
increment.  The  calculations  were  stopped  when  approached  zero, 

indicating  that  the  asymptotic  or  wet-bulb  temperature  has  been  approached. 
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APPENDIX  B 

CALCULATIONS  OF  PROPERTIES  IN  FILM 


The  computations  described  in  the  analysis  require  the  use  of  the 
thermal  conductivity,  viscosity,  density,  and  specific  heat  of  the  air- 
vapor  mixture  in  the  film  surrounding  the  droplet.  The  concentrations 
of  fuel  vapor  and  air  are  not  constant  across  the  film,  with  the  mixture 
being  the  richest  in  the  vicinity  of  the  droplet  and  leaner  away  from 
it.  To  a good  approximation  the  molal  ratio  of  fuel  vapor  to  air  equals 
the  ratio  of  the  partial  pressure  due  to  the  fuel  vapor  to  the  partial 
pressure  due  to  the  air  in  the  film.  Thus  the  largest  fuel-vapor  - air 
ratio  occurs  at  the  surface  of  the  droplet  and  is  essentially  equal  to 
the  ratio  of  the  vapor  pressure  at  the  droplet  temperature  to  the  partial 
pressure  due  to  the  air,  or  PfL/PaL*  This  ratio  becomes  zero  at  the 

outer  edge  of  the  film,  provided  there  is  no  interaction  between  droplets. 

It  has  been  shown  for  plane  surfaces  (ref.  8)  that  in  the  case  of  a 
fully  permeable  film,  where  mass  transfer  takes  place  in  both  directions 
and  equimolal  diffusion  exists,  the  partial  pressure  gradients  would  be 
linear  throughout  the  film.  The  temperature  in  the  film  which  in  most 
of  the  cases  varies  from  a maximum  of  Tf}>  the  air  temperature  at  the 

outer  edge  of  the  film,  to  T^,  the  droplet  temperature  at  the  inner 
edge,  would  follow  a linear  gradient  only  in  the  case  of  no  mass  trans- 
fer. In  the  case  under  consideration,  that  is,  a semipermeable  film 
with  mass  transfer  in  one  direction  only,  where  part  of  the  incoming 
heat  is  picked  up  by  the  diffusing  vapor,  the  linearity  of  both  the 
partial  pressure  and  temperature  gradients  in  the  film  ceases  to  exist. 

An  exact  evaluation  of  the  properties  in  the  film  surrounding  a 
spherical  droplet  would  render  the  equations  almost  impossible  to  inte- 
grate, and  furthermore  the  precision  with  which  some  of  the  properties 
are  known  probably  does  not  justify  the  work  required.  Thus  an  averaging 
technique  in  the  film  was  employed.  The  properties  then  have  been  cal- 
culated at  the  logarithmic  average  temperature  in  the  film  Tm  and  for 

a fuel-vapor  - air  mixture  of  constant  concentration  equal  to  one-half 
the  concentration  at  the  droplet  surface.  The  following  equations  were 
used: 


T, 


■m 


l°Se 


Tb 

Tl 


(29) 
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“m  - (l  - ?rVa  + 1% 


2pi 


T/ 


2P, 


T 


(50) 


pfL\  . pfL 


^ M1  " + 5f  ** 


(31) 


Pm  “ 


12RT, 


m 


rT 


12RT, 


m 


PfL 
2Pt 


(32) 


CPm  = gaCPa  + gpCPf 


°Pa  + » — °Pf 


Mm 


Mm 


(33) 


In  these  relations  the  pure  properties  of  the  fuel  vapor  and  air 
have  been  taken  at  the  logarithmic  mean  temperature  in  the  film. 

Figure  22  shows  the  properties  in  the  film  versus  droplet  tempera- 
ture as  calculated  for  n-octane  and  air. 


1 
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APPENDIX  C 

JUSTIFICATION  FOR  ASSUMPTION  OF  HIGH  FUEL  CONDUCTIVITY 
AND  UNIFORM  DROPLET  TEMPERATURE 


As  can  be  deduced  from  equation  (21),  the  calculations  performed 
in  this  project  assumed  a large  thermal  conductivity  through  the  liquid 
droplet  and  a resulting  uniform  temperature  throughout  the  droplet  at 
all  times.  An  estimation  of  the  validity  of  the  above  assumptions  can 
be  made  from  the  charts  prepared  by  Gurney  and  Lurie  (ref.  9;  discussion 
given  in  ref.  10)  for  determining  the  relationship  between  temperature 
and  time  at  various  points  in  spheres  and  other  shapes.  To  adapt  this 
method  to  the  cases  of  interest  here,  the  resistance  ratio  kjihjr0  was 

used,  where  kL  is  the  thermal  conductivity  of  the  liquid  fuel  evaluated 
at  the  temperature  of  the  droplet;  h^,  a coefficient  of  heat  transfer  to 

the  liquid  droplet  based  on  the  portion  of  heat  transferred  to  the  droplet 
that  appears  as  sensible  heat  (Ql)»  and  rQ,  the  radius  of  the  droplet. 

The  use  of  this  parameter  together  with  the  charts  permits  an  estimation 
of  the  possible  difference  between  the  temperature  at  the  surface  of  the 
droplet  and  the  temperature  at  its  center,  respectively.  The  Gumey- 
Lurie  charts  specify  that  the  temperatures  at  the  midpoint  and  at  the 
surface  would  differ  considerably  when  the  resistance  ratio  is  small  and 
only  slightly  when  it  is  large.  Referring  specifically  to  the  Gurney- 
Lurie  charts  for  spheres,  if  the  resistance  ratio  is  6 or  more,  the  tem- 
perature difference  is  negligibly  small;  while,  if  the  resistance  ratio 
is  1 or  2,  the  temperature  difference  is  appreciable  and  should  be  taken 
into  account.  However,  the  Gurney-Lurie  charts  were  based  on  the  assump- 
tions of  constancy  in  the  temperature  surrounding  the  surface  of  the  body, 
the  resistance  ratio,  and  the  thermal  diffusivity  in  the  body  under  con- 
sideration; and  It  should  be  realized  that,  under  the  conditions  of  an 
evaporating  droplet,  the  charts  should  be  used  for  estimation  purposes 
only. 


Values  for  the  resistance  ratio  are  plotted  in  figure  23  for  n-octane, 
benzene,  and  n-pentane  versus  heating-up  time  together  with  the  droplet 
temperatures  calculated  by  the  methods  previously  described  in  this  report 
for  an  air  pressure  of  1 atmosphere  and  an  air  temperature  of  500°  F. 
Figures  24(a)  and  24(b)  show  the  resistance  ratio  plotted  for  benzene  at 
a constant  air  pressure  of  1 atmosphere  and  a constant  air  temperature 
of  500°  F,  respectively. 

The  data  of  figures  23  and  24  indicate  that  the  assumption  of  high 
thermal  conductivity  in  the  liquid  droplet  is  unquestionably  good  for 
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high-volatility  fuels,  is  of  reasonable  validity  for  medium- volatility 
fuels,  and  should  be  further  investigated  for  the  low-volatility  fuels. 
If  internal  circulation  (ref.  7)  exists  in  t1  ; droplets  under  the  con- 
ditions of  interest  here,  it  could  effectively  increase  the  conductivity 
of  the  liquid. 
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Figure  4.-  Schematic  diagram  of  vaporization  equipment 


Figure  5»-  Vaporization- time  history  for  droplet  of  n-octane  with  initial 
radius  of  50  microns.  Ti  = 200°  F;  Ut  = 400  fps;  Pm  = 1 atmosphere; 
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Figure  8.-  Vaporization- time  history  for  droplet  of  n-octane  with  initial 
radius  of  100  microns.  T-i  = 50°  F:  T-n  = 100°  F. 
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(b)  At  constant  air  pressure  Pm  of  1 atmosphere 


= 500°  F 
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Figure  10.-  Temperature-time  histories  of  benzene. 

U-,  = 100  fps. 


(c)  With  constant  initial  radius  rol  of  50  microns 


Figure  11.-  Temperature- time  histories  of  n-pentane 

Ui  = 100  fps. 


oOOO 
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(c)  With  constant  initial  radius  rn-i  of  50  microns 


Figure  12.-  Tenrperature-time  histories  of  three  fuels. 

= 100  fpsj  rol  = 50  microns. 


Figure  12.-  Concluded. 
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(a)  Pip  = 0.5  atmosphere;  Tg  = 100°  P. 


Figure  13.-  Penetration-time  and  mass-time  histories  during  unsteady 
state  for  n-octane  at  nine  different  ambient  air  conditions. 
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(c)  PT  = 0.5  atmosphere;  TB  = 1,000°  F. 


Figure  13.-  Continued 
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(d)  Pip  = 1 atmosphere;  TB  = 100°  F. 


Figure  13,-  Continued 
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(e)  Pip  s 1 atmosphere;  Tg  = 5OO0  F. 


Figure  13.-  Continued 
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(f)  Piji  = 1 atmosphere;  TB  = 1,000°  F. 


Figure  13.-  Continued 
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(g)  P,j  = 5 atmospheres;  T-g  = 100°  F. 
Figure  15.-  Continued. 
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(i)  Pip  = 5 atmospheres;  T-g  = 1,000°  F. 
Figure  13.-  Concluded. 
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0 200  400  600  800  1,000  1,200 

Air  temperature,  TB , °F 


Figure  14 Comparison  of  calculated  and  experimental  wet-bulb  temperatures 

at  atmospheric  pressure. 


Temperature 


Figure  lp.-  Temperature-time  history  of  a droplet  of  n-octane. 


2  Room  air  temperature 

2 Droplet  placed  on  thermocouple 

3 Rapid  self-cooling  to  approximately  wet-bulb  temperature 

at  room  air  temperature 
k Gradual  approach  to  wet-bulb  temperature 

4 Hot  air  suddenly  swung  under  droplet 

5 Heating-up  or  unsteady- state  period  to  wet-bulb  temperature 

corresponding  to  hot  air 

6 Steady  state  at  wet-bulb  temperature 

7 Heating  up  of  thermocouple  to  air  temperature 


Figure  17.-  Effect  of  fuel  composition  on  heating-up  time  at  constant 

air  temperature. 
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Figure  18.-  Temperature,  velocity,  and  rate  of  mass  transfer  histories 
of  two  droplets  of  n-octane  of  50-micron  initial  radius.  Tp  = 50°  F; 

Up  = 100  fps;  Pip  = 1 atmosphere;  Tg  = 1,000°  F. 
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Figure  22  — Average  properties  of  air- vapor  mixture  in  film  at 
different  droplet  temperatures.  Air  pressure,  1 atmosphere 
air  temperature,  200°  F. 
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Figure  25.-  Resistance  ratio  and  droplet  temperature  for  three  fuels  at 
1-atmosphere  air  pressure  and  500°  F air  temperature . U]_  = 100  fps; 

roi  = 50  microns. 
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(a)  At  constant  air  pressure  of  1 atmosphere. 

Figure  24.-  Resistance  ratio  and  droplet  temperature  for  benzene 
Ui  = 100  fps;  rol  = 50  microns. 
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